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Using first-principles total-energy methods, we investigate Jahn- Teller distortions in III-V dilute 
magnetic semiconductors, GaAs:Mn and GaN:Mn in the cubic zinc blende structure. The results 
for an isolated Mn impurity on a Ga site show that there is no appreciable effect in GaAs, whereas, 
in GaN there is a Jahn- Teller effect in which the symmetry around the impurity changes from Td to 
T)2d or to C2t) . The large effect in GaN occurs because of the localized d* character, which is further 
enhanced by the distortion. The lower symmetry should be detectable experimentally in cubic GaN 
with low Mn concentration, and should be affected by charge compensation (reductions of holes 
and conversion of Mn ions to d^ with no Jahn- Teller effect). Jahn- Teller effect is greatly reduced 
because the symmetry at each Mn site is lowered due to the Mn-Mn interaction. The tendency 
toward ferromagnetism is found to be stronger in GaN:Mn than in GaAs:Mn and to be only slightly 
reduced by charge compensation. 

PACS numbers: PACS numbers: 75.50.Pp, 76.50.-fg, 72.25.Dc 



I. Introduction 

Electronics based upon the spin of the electron (spin- 
tronics) seeks to exploit the spin of charge carriers in 
semiconductor o^i^i^ . It is widely expected that new 
functionalities for electronics and photonics can be 
derived if the injection, transfer and detection of carrier 
spin can be controlled above room temperature in 
these dilute magnetic semiconductors (DMS)i.. Most 
of the work in the past has focused on InAsiMnSiSiS, 
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GerMnSLS^. Novel control of magnetism has already 
been achieved in these host materials. However, the 
reported Curie temperatures^ are too low to have signif- 
icant practical impact. Recently, there has been interest 
in wide bandgap semiconductors, such as GaN2&, which 
may exhibit higher Curie temperatureo^^i'^^i^^i'^^i^'^i'^^i^^. 
Since there has been tremendous progress on the growth 
of high-quality (Ga, Mn)N epitaxial layers^Si^l, GaN:Mn 
is a promising high Tc ferromagnetic semiconductor. 
Using the Zener model of ferromagnetism^^, Dietl et 
alt22. predicted that cubic GaN doped with 5 at.% 
of Mn and containing a high concentration of holes 
(3.5x 10^*'cTO~'^) should exhibit a Curie temperature 
exceeding room temperature. However, the mechanisms 
of ferromagnetism in DMS materials is still an open 
question. 

Theoretically, there are many proposals for the elec- 
tronic configuration which have focussed the ferromag- 
netic mechanism upon the coupling between the host 
p and Mn 3d state3^2iiSiiiii24^; suggestions include 
Mn 3d^-fhole induced ferromagnetism in GaAs:Mri2&, 
whereas in GaN:Mn the spin-spin interaction has been 
proposed to be driven by a double exchange mecha- 
nism involving of d-electrons in Mn Sd** states^. It has 
been suggested that 3d^, 3d^-|-holc, and 3d'* coexist in 
GaAs:Mn based upon evidence in the dilute regime that 
there is a Jahn- Teller— distortion associated with the 
3d^ stated. However, to our knowledge only one work 



has proposed that there should be a strong Jahn- Teller 
effect in the GaN:Mn system'*'^, and there have been no 
quantitative theoretical studies of the Jahn- Teller effect 
in either GaAsiMn or GaNiMn. In this paper, we report 
first principles calculations of the magnitude of the Jahn- 
Teller effect and the consequences for ferromagnetism in 
GaAs:Mn and GaN:Mn. 

The Jahn- Teller effect was first proposed to occur in 
open-shell molecules^, and there are many examples in 
impurity states in II- VI^ and III-\4l& semiconductors. 
The effect is caused by a distortion that lowers the sym- 
metry and leads to a splitting of a degenerate state that 
is linear in the magnitude of the distortion (Fig.l (a)). If 
the state is partially occupied, the total energy is always 
lowered by some distortion since all other contributions 
to the energy are quadratic in the distortion (Fig. 1(b)), 
and the total energy is minimum in a distorted configu- 
ration (Fig.l (c)). It may also happen that a large dis- 
tortion occurs leading to a new bonding configuration 
separated by an energy barrier, such as the DX center— 
and AX center—. Although the existence of a Jahn- Teller 
effect is determined only by symmetry, we have to do a 
quantitative calculation of the magnitude. If the splitting 
is less « 0.01 eV, we estimate that the effect is smaller 
than the effect of quantum fluctuation or temperature 
and therefore can be ignored. 

Our calculations were performed using the density- 
functional theory within the generalized gradient approx- 
imation (GGA) of the Perdew-Wang 91 form^. We 
used the Vanderbilt ultrasoft pseudopotentialsS and the 
Vosko-Wilk-Nusair interpolation for the correlation func- 
tional in the spin-polarized calculations, as implemented 
by the plane- wave total energy VASP code'^^. The cal- 
culated lattice constants are a = 5.751 A for GaAs and 
a = 4.542 A for GaN, and all calculations with super- 
cells are done keeping the supercell lattice vectors fixed 
as multiples of the primitive lattice vectors. All cell- 
internal structural parameters are fully relaxed until the 
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forces are converged to within 0.05 eV/A. The cutoff en- 
ergy for the planewave expansion is 170 eV for calcula- 
tions involving Mn in GaAs and 270 eV for Mn in GaN, 
with check using 400 eV. For self-consistent total energy 
calculations, we used a 64- atom supercell and a 4x4x4 
Monkhorst-Pack fc-point mesh (which has been verified 
to be sufRcieniiSS.) , and selected cases were checked with 
a 6x6x6 fc-point mesh. Density of states plots were made 
using a finer 8x8x8 fc-point mesh. 

I. ISOLATED MN IMPURITES 

In order to study isolated Mn impurites, we have car- 
ried out calculations on 64 atom cells with one Mn sub- 
stituted for a Ga atom in GaAs and in GaN. We have 
calculated the total energy and the eigenvalues of the 
Kohn-Sham hamiltonian for various cases. First we have 
considered the ideal geometry with all atoms at the ideal 
positions of the GaAs or GaN tetrahedral lattices, and 
breathing distortions in which the bond lengths change 
but the symmetry is constrained to remain T^. We have 
also considered three different lower symmetry distor- 
tions as shown in Fig. 2: (a) to T)2d (b) to C2v, 
(c) Td to Csi,. All cases can be compared if we define 
symmetry-adapted variables. In order to disentangle the 
effects of the distortions, we define the positions of the 
four nearest-neighbor N atoms around a Mn atoms to be 
R?, z = 1,4 in the ideal tetrahedral structure, and the 
displacements of the neighbors relative to the Mn, as: 

AR, = (Ax,,Ay,,Az,) (1) 

The displacement of the neighbors can be projected into 
a two parts. One part is a symmetric radial "breathing" 
component, which preserves the tetrahedral symmetry, 
given by 

ARf = AR^ R°, (2) 

where the magnitude AR^ is easily extracted for a given 
displacement pattern AR^ using 

Ai?^ = i^R°.AR,. (3) 

i 

The remaining parts of the displacement are symmetry- 
breaking Jahn- Teller distortions, given by 

ARf^ = AR, -ARf. (4) 

In one set of calculations, we have considered only 
breathing, varying the magnitude of ARf and forcing the 
symmetry to remain Td. Comparing to ideal tetrahedral 
positions, for GaN all four N nearest neighbors of the Mn 
atom move closer to the Mn by the amount ARf = 0.034 
A and the energy decreases by AE = 0.05 eV. For GaAs 
all four As nearest neighbors of the Mn atom move closer 
to the Mn by only a small amount ARf = 0.007 A and 



the energy decreases by AE = 0.008 eV. The difference 
between GaAs and GaN is readily explained by the dif- 
ferent size of As and N atoms. 

In the remaining calculations the atoms are allowed to 
distort and lower the symmetry. In the case of GaAs 
the energy was never found to decrease. Therefore, we 
conclude that there is no appreciable Jahn- Teller effect in 
GaAs. However, in GaN the energy is found to decrease 
substantially. The decrease in the total energy for the 
three different symmetries is listed in Table I. There we 
have defined the energy relative to the minimum energy 
of Td symmetry as described above. In the distortion 
shown in Fig. 2 (a), the four nearest neighbors of Mn 
atom move to lower the symmetry from T^ to T)2d with 
displacements (|Axi| — \Ayi\ ^ lAz^l). The distortion 
results in a 0.10 eV lower energy in GaN:Mn (see Table 
I the T)2d column). (This energy difference changes by 
only 6 meV if the energy cutoff is increased from 270 eV 
to 400 eV). In the Fig. 2 (b), two neighbors of Mn atom 
move along [110] and another two neighbors of Mn atom 
move along [110] direction. This atomic configuration 
with 0.08 eV total energy lowering in GaN:Mn (Table I, 
C2v column). However, there is not an obvious effect for 
Fig. 2 (c) distortion, in which Mn atom move along [111] 
direction, as shown in the column for Csy in Table I. 

From these calculations, we conclude that a strong 
Jahn- Teller effect should be observed for low concentra- 
tions in GaN (for and C2v) but not in GaAs:Mn. 
Furthermore in GaN:Mn the distortion is slightly favored 
for D2d symmetry; however, the energies are close enough 
that either distortion may occur in actual systems. We 
note that the lowering of energy from breathing relax- 
ation alone (0.05 eV) is smaller than the lowering of the 
energy (0.10 eV and 0.08 eV relative to breathing relax- 
ation) for the Jahn- Teller distortion (both D2d and C2t,). 
This shows clearly the importance of the Jahn- Teller ef- 
fect. In addition we have analyzed the magnitudes of the 
components of the T)2d displacement for the lowest en- 
ergy. The magnitude of the breathing component AR^ 
is 0.020 A and the magnitude of the Jahn- Teller displace- 
ment is |AR/-^| — 0.068 A for each of the four neighbors. 
Thus the magnitude of the components of the displace- 
ment, also show that the Jahn- Teller effect, is larger than 
the breathing for Mn in GaN. 

We interpret as a Jahn- Teller distortion caused by the 
partial occupation of the t2d state which is split due to 
the lower Z?2d symmetry. The Mn triplet t2d of 3d energy 
level splits into singlet and doublet states at F point. The 
splitting increases approximately linearly with the mag- 
nitude of the atomic displacements, as expected. At the 
minimum energy position, the splitting energy AEt^^ is 
0.23 eV. We have studied the stability of the distorted 
state by varying the magnitude of the distortion in vari- 
ous ways. 

We also considered large distortions in Fig. 2 (b) and 
(c) that might lead to qualitative rebonding of the atoms, 
e.g, two N (or As) atoms forming nearest neighbors as in 
an AX center (like Fig. 2 (b)) or large displacements of 



3 



the Mn atom along the [111] direction as in a DX center 
(Hke Fig. 2 (c)). However, no configuration was found to 
be stable or even metastable state in either GaAs:Mn and 
GaN:Mn. 

Figure 3 shows the partial density of states (DOS) of 
the Mn 3d states in GaAs:Mn (Fig. 3 (a)) and cubic 
GaN:Mn in symmetry (Fig. 3 (b)) and the Jahn- 
Teller distortion of D2d symmetry (Fig. 3 (c)). The re- 
sults for the undistorted cases (Fig. 3 (a) and 3 (b)) are 
very similar to those found by other groupsSiSi^SiS. In 
GaAs:Mn, we find that Mn 3d states is 2 eV below host 
valence band minimum (VBM) and the t2d level is lower 
than the e^ level (see Fig. 3 (a)). On the other hand, 
in GaN:Mn, the Mn 3d state is above host VBM and 
the t2d level is higher than (see Fig. 3 (b)). Because 
of this difference between GaAs:Mn and GaN:Mn, there 
are different consequences for the Mn 3d states when the 
symmetry changed from to T)2d- the t2d energy level 
has no splitting (< 0.01 eV) in GaAs:Mn, but a large 
splitting (0.23 eV maximum at F point) in GaN:Mn, with 
the Fermi level in the gap as shown in Fig. 3 (c). 

For the case of an isolated Mn substituted for a Ga 
atom in GaN, we also considered the charged state with 
an added electron. This charge state is expected when 
the system is compensated, reducing the number of holes 
and leading to Mn atoms in the 3d^ state. Since this 
is a symmetric closed-shell state, the Jahn- Teller distor- 
tion should disappear. Indeed, our calculations bear this 
out. We carried out calculations in which the atoms are 
constrained to have T^ symmetry (pure breathing relax- 
ation) and other calculations in which the atom positions 
are relaxed starting with distorted initial positions having 
D2d symmetry. In all cases the total energy is the same 
to within 0.02 eV, the splitting of t2d is very small (about 
0.01 eV), and the atoms relax to positions near T^ sym- 
metry . In principle, all calculation should have the same 
total energy. We consider the differences to be negligible 
and these results can be considered as numerical tests 
showing the accuracy of our calculations. Clearly, this 
supports our conclusion that an isolated neutral Mn sub- 
stituted for Ga has a Jahn- Teller distorted ground state 
in GaN, and that compensation can decrease the effect 
causing some Mn to in the 3d^ state with no Jahn- Teller 
distortion. 



II. INTERACTIONS OF MN-MN PAIRS 

Although we have shown that isolated Mn impurities 
undergo a large Jahn- Teller distortion in GaN, further 
studies are required to establish the effects, if any, upon 
the properties of GaN:Mn alloys in interesting concen- 
tration ranges for magnetic semiconductors. Since ferro- 
magnetism is due to interactions between the Mn atoms, 
we must consider the effect upon Mn-Mn interactions for 
Mn atoms separated by typical distances found in the al- 
loys. The Jahn- Teller effect would be expected to cause 
a reduction in the tendency toward ferromagnetism since 



it leads to a splitting of the states. The splitting will 
reduce the interactions between Mn pairs since hoping 
requires an extra energy cost. On the other hand, the 
interaction between the Mn atoms may be so large that 
it dominates over the Jahn- Teller effect. 

In order to study Mn-Mn interactions, we have car- 
ried out calculations on 64 atom cells with two substitu- 
tional Mn atoms at various distances and in different spin 
states. We selected 2 Mn atoms separated by distances 
\/2a and V2a/2, where a is the cubic GaAs or GaN lattice 
constant. For GaAs case, for separation v^a, the total 
energy difference between ferro- and antiferro-magnetic 
spin states AEaf—0-'2'2 eV/Mn-pair, which is compa- 
rable with 0.2 eV/Mn-pair in Ref.'^?, and 114 meV/Mn 
(0.228 eV/Mn-pair) in Ref.^^. For the GaN case, we have 
considered different starting configurations of the atoms, 
in one set of calculations starting with all atoms in the 
ideal zinc blende positions and in a second set starting 
with Jahn- Teller distorted states (the atoms around the 
Mn-Mn pairs are placed in distorted configurations). The 
final positions of the atoms after relaxation are the same 
for the two cases, showing that checked that final con- 
figuration was independent of the starting point. For 
separation \/2a/2, the energy difference is AEaf—0.36 
eV/Mn-pair, this value is good in comparing to similar 
workSi^, which is 188 meV/Mn and 156 meV/Mn. For 
separation \^a, the AEaf=0-50 eV/Mn-pair— , while a 
similar workS gave 161 meV/Mn. 

The energy levels for the Mn 3d states in the gap show 
that the interactions between the two Mn atoms in a pair 
is indeed larger than the splitting caused by the Jahn- 
Teller effect on the individual Mn atoms. In the ferro- 
magnetic state the d energy levels are split by the Mn- 
Mn interactions by amounts that are much larger than 
the splitting due to the Jahn- Teller effect. The maximum 
width of the d-bands in our supercell calculation is 0.4 
eV for the Jahn- Teller splitting as shown in Fig. 3 (c). In 
contrast, the maximum width of the d-bands is 1.4 eV for 
a ferromagnetic Mn-Mn pair separated by -\/2a/2 with all 
atoms relaxed. Furthermore, the Mn-Mn interactions are 
present even if all atoms are in ideal positions, in which 
case the maximum width is only slightly changed to 1.3 
eV. Similar results are found for Mn-Mn pairs separated 
by \/2a.. These results show that Jahn- Teller distortions 
do not have a large affect upon the magnetic interactions 
between Mn atoms at distances expected in actual ferro- 
magnetic semiconductors. 

It is interesting also to consider the realistic case with 
charge compensation, in which some of the t2d levels are 
filled. If the states are localized near the Mn atom as 
in GaN:Mn, each added electron can be interpreted as a 
conversion of a d^ into a d^ state. Since a d^ state is a 
spatial singlet with no degeneracy except spin, no Jahn- 
Teller effect will occur. We also calculated -1 charged 
states for Mn pair separated with a and \/2a/2 GaN 
lattice constant. All the results as well as neutral Mn-Mn 
pair case are listed in Table II. We found that the total 
energy difference between ferro and anti-ferro magnetic 
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states are decreased in charged states, such as for ^/2a/2^ 
AEaf changes from 0.36 eV/Mn-pair to 0.30 eV/Mn- 
pair , and for y/2a case, AEaf changes from 0.30 eV/Mn- 
pair to 0.22 eV/Mn-pair (sec Table II). This result show 
that charge compensation will decrease the tendency for 
ferromagnetism, as expected. 

III. CONCLUSIONS 

In summary, we have studied the Jahn- Teller distor- 
tion in GaAs:Mn and GaN:Mn. Our rcsidts show that 
a strong Jahn- Teller distortion should happen in uncom- 
pensated GaN:Mn at low concentrations where the Mn 
impurities arc isolated. The lowering of the energy is 
due to the splitting of the t2d states of the localized d- 
electrons on the Mn 3d^ ion, leading to an energy gap. 
There are two possible symmetries C2v and D2d, with the 
latter having the lowest energy. In the presence of charge 
compensation, the Mn d states are filled leading to filled 
shell spherically-symmetric 3d^ ions and the Jahn- Teller 
effect disappears. These effects should be observable ex- 
perimentally. 

In contrast, in GaAs:Mn the Mn 3d states are primarily 
3d^ with a hole in the GaAs valence band. This state is 
only weakly coupled to the distortions and the tendency 
for a Jahn- Teller distortion is a negligible effect. 

In order to study the effects upon magnetism, we car- 
ried out calculations on Mn pairs. In agreement with 



other work, we find Mn-Mn interactions to lead to fer- 
romagnetism in both GaAs and GaN, with larger inter- 
action in GaN. The interaction between Mn-Mn pairs at 
realistic distances is sufficiently large that it dominates 
over the Jahn- Teller effect. The interactions between Mn 
atoms is not greatly affected by lattice relaxations and 
there is always a clear tendency for ferromagnetic align- 
ment of Mn pairs. In realistic cases with charge compen- 
sation, ferromagnetism is still favored: even with 100% 
compensation the ferromagnetic interactions are reduced 
by only 20%. Finally, even though our calculations are 
for the cubic structure, the conclusions on magnetic inter- 
actions should carry over to the wurtzite structure since 
they do not depend upon detailed positions of the atoms 
and the ferromagnetic state persists whether or not there 
are distortions. Thus our results support the conclusion 
that GaN:Mn holds promise as a ferromagnetic semicon- 
ductor. 
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FIG. 3: Partial density of states and schematic energy level 
diagrams of Mn 3d states in (a) GaAs:Mn, (b) GaN:Mn in Td 
symmetry without a Jahn- Teller distortion (no energy level 
splitting), (c) With Jahn- Teller distortion of D2d symmetry, 
which shows the splitting of the t2d d-states. 
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TABLE I: Comparison of total energy among four different 
symmetries in GaN:Mn. Here the total energy Etot of 
symmetry (with breathing relaxation) is set to zero. 



Symmetry 




D2d 


C2v 




Etot (eV) 


0.00 


-0.10 


-0.08 


-0.02 



TABLE II: Total energy difference AE between ferro- 
(FM) and antiferro- (AFM) magnetic spin states for a Mn-Mn 
pair in neutral and -1 charged GaN:Mu. Here AE^^ =EfJt'^' - 
Etot' 1 •\/2o/2 and \/2a are two separations of the Mn-Mn pair, 
and a is lattice constant of cubic GaN. 



System 



Neutral GaN:Mn (-1) charged GaNiMn 



4 



Spin Up 



2 







2 



4 




Partial DOS projecl 



Spin Down 

}26 ^ - 



on Mn-3d states 




Partial DOS projected on Mn-3d states 



4 



Spin Up 




Partial DOS 



Spin Down 

^2d ^ 



on Mn-3d states 



